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O ■ ABSTRACT 

The Transit Timing Variations (TTVs) can be used as a diagnostic of grav- 
itational interactions between planets in a multi-planet system. Many Kepler 
>• . Objects of Interest (KOIs) exhibit significant TTVs, but KOI-142. 01 stands out 

among them with an unrivaled, ~12-hour TTV amplitude. Here we report a 



C^I , thorough analysis of KOI-142. Ol's transits. We discover periodic Transit Du- 

• . ration Variations (TDVs) of KOI-142. 01 that are nearly in phase with the ob- 

O ■ served TTVs. We show that KOI-142. Ol's TTVs and TDVs uniquely detect a 

non-transiting companion with a mass ~0.7 that of Jupiter (KOI- 142c). KOI- 
142. Ol's mass inferred from the transit variations is consistent with the measured 
transit depth, suggesting a Neptune class planet (KOI- 142b). The orbital period 



;_i . ratio Pc/ Ph = 2.03 indicates that the two planets are just wide of the 2:1 reso- 

- -' nance. The present dynamics of this system, characterized here in detail, can be 

used to test various formation theories that have been proposed to explain the 

near-resonant pairs of exoplanets. 



1. Introduction 

The two methods used so far to discover the majority of exoplanets, transit observations 
(TOs) and radial velocity (RV) measurements, have known limitations. The best current 
RV sensitivity of ~1 m s~^ allows us to detect planets down to ~1 Neptune mass at 1 AU, 
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assuming a bright, solar-mass host star and sin J ~ 1, where / is the inchnation of planet's 
orbit relative to the sky plane. The TO method, on the other hand, can be used to detect 
smaller planets (Borucki et al. 2011), but requires that sin J ~ 1. Most multi-planet systems 
discovered from TOs are therefore implied to be coplanar to within ^1° (Lissauer et al. 
2011a). 

Ideally, we would like to use the detection statistics of planets and planetary systems 
from TOs and RVs to build a population model that globally describes the distribution 
of planet properties (such as the planet occurrence depending on the host star mass and 
metalicity, planetary mass and size distributions, multiplicity, distributions of orbital spacing, 
eccentricity and inclination, etc.) in the Galaxy. Much of this information, however, is 
difficult to obtain from TOs and RV measurements alone. The problem is rooted in the fact 
that these methods are blind to certain classes of planets and planetary systems, and often 
allow us to obtain only approximate properties of the detected systems. 

The Transit Timing Variation (TTV) method can help resolve some of these issues. The 
TTVs occur when gravitational perturbations produce deviations from strictly Keplerian 
orbits, such that the spacing of individual transits is not exactly periodic (Miralda-Escude 
2002, Agol et al. 2005, Holman & Murray 2005). The TTVs have been used to confirm some 
of the transiting planet candidates from Kepler (e.g., Holman et al. 2010, Lissauer et al. 
2011b), detect and characterize non-transiting planets (Nesvorny et al. 2012), and search for 
moons (Kipping et al. 2012, 2013). The TTV analysis often provides specific information 
about the system that is unavailable from TOs or RVs alone. 

Here we discuss the Kepler Object of Interest (KOI) 142.01. We show that the observed 
TTVs can uniquely be fit by a sub- Jovian mass planet near the exterior 2:1 orbital resonance 
with the transiting Neptune-class planet. The Transit Duration Variations (TDVs) produced 
by the interaction of two planets are also detected. We use the TTVs and TDVs to investigate 
the near-resonant dynamics of the KOI-142 system. The methods and results are described 
in Sections 2 and 3. The broader implications of our work are discussed in Section 4. Efforts 
such as these can can help us to understand the formation of the near-resonant planetary 
systems. They can also provide important guidance in our pursuit of the planet population 
model. 



2. Method 

The TTVs of KOI-142.01 were identified by Ford et al. (2011, 2012) and Steffen et al. 
(2012). Mazeh et al. (2013) updated KOI-142.01 's TTVs using the first twelve quarters of 
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the Kepler data, highlighted the large TTV amplitude (~12 hours), and suggested that this 
'king of TTVs' can host one or more additional planets. 

We downloaded the publicly available data for KOI-142 from the Mikulski Archive for 
Space Telescopes (MAST), which included short-cadence (SC) data from quarters 5-14 and 
long-cadence from quarters 0-14. Throughout our analysis, we make use of the Photometric 
Analysis (PA) time series. Using a polynomial ephemeris fitted to the Mazeh et al. (2013) 
transit times, we extracted each transit epoch with ±0.5Pp worth of data either side (Fig. [T]) 
and applied the CoFiAM (Cosine Filtering with Autocorrelation Minimization) detrending 
algorithm, which is described in detail by Kipping et al. (2013). 

Although we direct those interested to Kipping et al. (2013) for details, we briefly de- 
scribe CoFiAM. The algorithm is essentially a high-pass, low-cut filter which removes periodic 
components with a timescale greater than that of a pre-defined protected timescale. For this 
analysis, the protected timescale was chosen to be that of three times the transit duration to 
ensure the transit shape is minimally distorted by our detrending. CoFiAM is applied to each 
transit separately, after first removing discontinuous features by eye and applying a mov- 
ing median outlier filter. For each transit, CoFiAM explores typically dozens of permissible 
detrendings by permuting the maximum allowed harmonic order from the maximum (corre- 
sponding to the protected timescale) to the minimum (corresponding to the entire baseline 
of the time series). CoFiAM then selects the order which minimizes the autocorrelation on a 
30 minute timescale, as determined by the Durbin- Watson metric, d. 

After detrending the data, the next step is to fit a transit model to the cleaned, normal- 
ized photometry. For this purpose we use the Mandel & Agol (2002) algorithm to model the 
transits assuming a circular orbit and the multimodal nested sampling algorithm Multi- 
Nest (Feroz et al. 2009, 2011) for the regression. Due to the presence of large dynamical 
variations, we wished to investigate the possibility of transit duration variations (TDVs) in 
addition to TTVs. This requires allowing all of the basic transit parameters to vary for each 
event and so we regress every epoch individually. In total, we detrended and fitted 105 tran- 
sit epochs although 5 of these did not converge due to data gaps leading to partial transits. 
The M = 100 well-fitted transits lead to M joint-posteriors for the transit parameters. Our 
"best-fit" transit times and durations are computed by marginalizing each epoch's time of 
transit minimum (r) and duration defined as the time it takes for the planet's center to cross 
the stellar disc to exiting under the same condition (T, see Kipping 2010 for further details 
on this parameter). The final values are the median of each distribution and the ±34.1% 
quant iles and are shown in Figure ID 

The M posteriors may also be used to derive estimates of the mean transit parameters, 
such as Rp/R^:, a/R^, etc. For example, we consider a posterior of the mean Rp/R^ made of 
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N realizations, labeled i = 1 ^ N. For the i^^ realization of the mean Rp/R^, we calculate 
this as simply the sum of the i^^ realization of each of the M transit epochs divided by M. 
N is limited by the smallest length of the M joint posteriors, which was N = 14034. We 
provide the median and associated ±34.1% quantiles of these mean transit parameters in 
Table [H 

Next, we investigate dynamical solutions for the large TTVs and TDVs and in what 
follows we utilize the best-fit TTVs and TDVs and their associated uncertainties. We tested 
whether the measured TTVs and TDVs are consistent with gravitational perturbations from 
a planetary or stellar companion of KOI-142.01, and whether a unique set of parameters can 
be derived to describe the physical and orbital properties of that companion. We examined 
orbits with periods between 1 day and 10 years, including the cases of highly eccentric and/or 
retrograde orbits. 

The dynamical fits were obtained with a new code based on a symplectic iV-body inte- 
grator known as swif t_rmvs3 (Levison & Duncan 1994). The code computes the mid-transit 
times by interpolation. First, the transiting planet is forward propagated on the ideal Keple- 
rian orbit starting from the position and velocity recorded by swif t_rmvs3 at the beginning 
of A^-body time step. Second, the position and velocity at the end of the time step are 
propagated backward (again on the ideal Keplerian orbits). We then calculate a weighted 
mean of these two Keplerian trajectories such that progressively more (less) weight is given 
to the backward (forward) trajectory as the time approaches the end of the time step. 

The method described above is efficient^ The required transit timing precision was 
achieved by setting the timestep to ^1/20 of the inner orbit period. The TTVs were com- 
puted relative to a linear ephemeris. The dynamical fits to the TTVs were obtained by 
minimizing 

M 

XTTV = Y.^^to,,-Stc,jf/a] , (1) 

i=i 
where M = 100 is the number of transits, 6to.j and Stcj are the observed and calculated 
TTVs, and aj is the uncertainty of Stoj- 

The transit duration of each transit was determined from the impact parameter and the 
projected transit speed at mid transit. The dynamical fits to the TDVs were obtained by 



^The A^-body code is nearly as fast as the perturbation method (Nesvorny & Morbidelh 2008, Nesvorny 
& Beauge 2010) that was used to detect the non-transiting planetary companion of Kepler-46b (previously 
known as KOI-872b; Nesvorny et al. 2012). Unlike the perturbation method, however, the iV-body code can 
account for the resonant and near-resonant transit variations, which is important for KOI-142. 
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mininiizmg 

M 

xIdv = E(^o.-Tc,)VS?, (2) 

i=i 

where Toj and Tcj- are the observed and calculated TDVs, and Sj is the uncertainty of 
To J. The simultaneous dynamical fits to the TTV and TDVs were computed by minimizing 

X = Xttv + Xtdv 

We tested three dynamical models: (1) coplanar orbits and circular orbit of the transit- 
ing object (model C), (2) inclined orbit of the companion and circular orbit of the transiting 
object (model X), and (3) a general case where the transiting object has an eccentric orbit, 
and the companion orbit can be inclined and eccentric (model £). Model C has 6 parameters 
(two masses and four orbital parameters of the companion), model X has 8 parameters (six of 
model C plus the inclination and nodal longitude of companion's orbit), and model £ has 10 
parameters (eight of model X plus the eccentricity and pericenter longitude of the transiting 
planet). 

The remaining four parameters, namely the semimajor axis, inclination, nodal longi- 
tude and mean longitude of the transiting object at a given epoch, were held fixed. The 
semimajor axis was computed from the orbital period (obtained from a linear ephemeris). 
The inclination and nodal longitude were set from the impact parameter h. We used the 
transit reference system (Nesvorny et al. 2012), where the nodal longitude of the transiting 
object is 270°. The mean longitude was set so that the first transit occurred at the reference 
epoch (to = 2454954.62702 BJDutc)- 

To compute the errors on the parameters fitted in the dynamical fits, we exploit the 
posterior files derived for r and T for each transit epoch. Specifically, we draw the i^^ sample 
from the ensemble joint posterior for every epoch to create a fair realization of the TTV and 
TDV for this i^^ draw. We then repeat the process of the dynamical fits described above for 
the best-fit TTVs/TDVs. Since the draws do not technically have an error and the errors of 
the best-fit values were found to be approximately equal anyway, we use equal weighting in 
the subsequent regression. We then end up with a vector for the dynamical parameters best 
describing the z*^ draw of the joint posterior. This is repeated for 1000 draws (we did not do 
all 14034 available realizations due to computational constraints) allowing us to compile a 
joint posterior for the dynamical parameters (Figs. |2]and|3]). We found that 45 realizations 
resulted in unphysical results and these were discarded accordingly. The final posteriors were 
marginalized for each parameter of interest and we derived the tabulated values of Table [2] 
by quoting the median and the associated ±34.1% quantiles. 

Physical stellar parameters were derived by matching stellar evolution isochrones to the 
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observable stellar properties. Two spectra of KOI-142 were obtained using the Fibre-fed 
Echelle Spectrograph (FIES) at the 2.5 m Nordic Optical Telescope (NOT) at La Palma, 
Spain (Djupvik & Andersen 2010) on 14 and 15 July, 2011. We used the medium-resolution 
fiber with a resolving power of A/AA ~ 46,000 and an exposure time of 11 to 16 minutes 
yielding a SNR per resolution element of 21 to 28. We used SPG (Buchhave el al. 2012) 
to determine the stellar parameters of the host star, yielding an effective temperature of 
Teff = 5513±67K, a surface gravity of log c/ = 4.50±0.12, ametallicity of [m/H] = 0.37±0.08 
and a projected rotational velocity of vsini = 1.8 ± 0.5 km s~^. The relatively low SNR 
spectra are at the limit of what SPG requires to extract reliable stellar parameters, which 
may not be accurately reflected in the formal uncertainties but is revealed by the relatively 
small percentage of draws matching parameter sets which are allowed by stellar models (see 
below). 

Rather than use logg as a luminosity indicator, we opt to use p* derived from the transit 
light curve. Although one might suppose we have 14034 realizations of p*, we actually have 
14034 realizations of p*,circ since we assumed a circular orbit in the original transit fits. The 
true stellar density may be easily derived using the simple correction p* = p*,circ/^ where 
\l/ = (1 = e sin 07)^(1 — e^)~^/^ assuming e is not large (see Kipping 2011 for the derivation). 
Since we only have 955 dynamical realizations comprising the posteriors of e and u, we are 
only able to derive 955 realizations of ^, from which we derive "^ = 1.18321qqq^5. Using the 
corrected 955 fair realizations of p*, we draw a random normal variate for T^s and [Fe/H] 
determined by SPG and then match each of these 955 draws to the Yonsei-Yale theoretical 
stellar models (Yi et al. 2001). 

We find that the joint posterior of {p*,Tefr, [Fe/H]} lies on the edge of the permissible 
range allowed by the stellar models, with only 244 draws matching to a model. The draws 
that match tend to have lower metalicity, density and temperature than the median values 
from the full sample of 955 draws. From these trials, we marginalize over each parameter's 
posterior to estimate the physical parameters for the host star KOI-142. The final stellar 
parameters, defined as the median and their ±34.1% associated quantiles, are provided in 
Table [31 We use these to compute the physical parameters of planets b and c too, as shown 
in Table H 



3. Results 

KOI-142. 01 shows transits with a period of P ~ 10.95 days. The transit timing is 
modulated with a TTV period of Pttv — 630 days. The TTV amplitude exceeds 10 hours 
and appears to be changing with time (Fig. S^). We detect periodic TDVs nearly in phase 
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with the measured TTVs (Fig. |Dd). The TDV amphtude is much smaller than that of the 
TTVs. The large Pttv/P ratio (~58) suggests that the TTVs may be related to secular, 
resonant or near-resonant perturbations from a companion, rather than to the short-periodic 
effects. 

Our detailed dynamical modeling of the measured transit variations uniquely detects a 
non-transiting planetary companion near the exterior 2:1 resonance with KOI-142.01. The 
uniqueness of the fits stems from the successful modeling of both the large near-resonant 
TTVs shown in Fig. H^, and the short-periodic 'chopping' produced by the orbital conjunc- 
tions between planets (Fig. [5]). The retrograde orbits near the 2:1 resonance, for example, 
can be ruled out because they lead to a chopping pattern that is inconsistent with data 
(conjunctions occur too often). Below we discuss the results based on the general dynamical 
model S. 

Our best TTV-only fit gives Xttv — 65.7 for 90 degrees of freedom (DOF), while the 
best simultaneous fit to the TTVs and TDVs gives x^ = 182 for 190 DOF. The two fits 
are consistent with each other (within errors). They fit data very well (Fig. H]). All other 
solutions, including the highly inclined or retrograde orbits, can be ruled out because they 
give Xttv > 1000. Interestingly, the TDVs expected from the TTV-only fit are very similar 
to those obtained from the simultaneous fit (Fig. HJd). This shows a good consistency of the 
identified solution. The best fit parameters and their errors are listed in Tables [2] and HI 

We find that the companion mass is 6.3 x 10~^ M*, where M* = 1.02 Mq (Table 3). 
This suggests a planet with mass ~0.67 Mj, where Mj is the mass of Jupiter, or about 2.2 
Saturn masses. The KOI-142.01's mass inferred from the TTVs is < 5 x 10~^ M*, implying 
a sub-Neptune mass planet. This is consistent with the radius ratio i?/-R* = 0.039 inferred 
from the transit analysis. With R^, = 0.96 Rq, where Rq is the Sun's radius, this would give 
a planetary radius similar to that of Neptune's. We therefore confirm KOI-142 as a system 
of two planets, hereafter KOI-142b and KOI-142c. 

KOI-142b's and KOI-142c's orbital periods are Pb ~ 10.95 days and P^ ~ 22.34 days, 
respectively. The orbital period ratio of the two planets is therefore Pc/Pb = 2.04 (or ~ 2.03 
if averaged over semimajor axis oscillations; Fig. ^, indicating orbits just wide of the 2:1 
mean motion resonance. This orbital configuration is relatively common among the KOIs 
and confirmed Kepler planets (Fabrycky et al. 2012). Our long-term integrations of the 
system show that the orbits are stable on (at least) Gyr timescales. 

The two resonant angles, ab = 2Ac — \b — '^b and Uc = 2Ac — A^ — Wc, where A's and 
Clt's are the mean and pericenter longitudes, circulate in a retrograde sense with a period 
of P2:i — 630 yr. The associated eccentricity variations of KOI-142b's orbit produce the 



observed TTVs and TDVs. A detailed interpretation of TTVs can be obtained from Lithwick 
et al. (2012). Here we just point out that given the proximity to the 2:1 resonance, at and 
Gc are not simple linear functions of time. The transit variations are therefore not strictly 
sinusoidal, as noted by Mazeh et al. (2013). 

As for the TDVs, the transit duration can be approximated by T = A/V^, where A = 
2_R*Vl — &^ is the length of the transit chord, and V = nf,af,(l + Cb cos Wb) / ^/ 1 — el is the 
projected speed. Here, n^ = 271 /Pb denotes the orbital frequency of KOI-142b. We find that 
the main contribution to the TDVs of KOI-142b comes from the variation of 1/(1 + e cos cc), 
so that 

_,^^_2wrH?,,^_ (3) 

where kb = eb cos Wb and Skb denotes the variation of kb. The TDVs produced by the variation 
of Ubttb are much smaller. The vertical TDVs produced by the variation of Vl — 6^ become 
important on longer timescales (see below). 

Interestingly, despite their very different masses, KOI- 142b and KOI- 142c have similar 
orbital eccentricities (mean Cb = 0.064 and Cc = 0.055; Fig. [6]). This is an interesting 
constraint on theories that attempt to explain the near-resonant planet pairs (Fabrycky et 
al. 2012) by tidal migration (Terquem & Papaloizou 2007, Lithwick & Wu 2012, Batygin & 
Morbidelli 2012) or by planet-disk interactions (Baruteau & Papaloizou 2013). The relatively 
large eccentricity of KOI- 142c needs an explanation. 

The lines of apses of the two planets are offset so that Azu = vJc — vjb — 180° at the 
reference epoch, but Aro is not stationary, because Wb circulates in a prograde sense with 
the period P^ — 20 years, while the secular drift of Wc is much slower. The proximity of Ace 
to 180° at the reference epoch is therefore probably coincidental. The secular precession of 
Wb on the P^ period correlates with the ~20-yr modulation of the resonant oscillations of 
Cb seen in Fig. |6b. 

We find a bimodal distribution of Vt^ with Vt^ ~ 90° or Vie ~ 270°. The apparent absence 
of transits of KOI-142c implies that ic > 1.5°. Our best fits suggest that 2° < ic < 6°. While 
the two orbits are therefore approximately coplanar, the mutual inclination is significant, 
explaining why both planets' transits were not detected. The nodal precession period is 
Pq c^ 100 yr. The secular oscillations of ib happen on this timescale (Fig. [Hb)- The 
variations of ic have much smaller amplitude because Mc/Mb > 13. 

The secular dynamics of KOI- 142b has interesting implications for the future behavior 
of TTVs and TDVs. As shown in Fig. [7^, the TTV signal is expected to be modulated on 
the Pro period. Moreover, the secular changes of ib and Qb will profoundly affect the transit 
duration (Fig. [Tjo) and the impact parameter (Fig. Wp)- Depending on whether Qc ~ 90° 
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or Qc ~ 270° at the present, the transit impact parameter bb is expected to decrease (for 
Q^ ^ 90°) or increase (for Qc — 270°) over the next few years. Monitoring bh over the next 
decade will thus help to remove the Qc degeneracy, and identify the correct value of Qc- 

In any case, we predict that KOI-142b should stop transiting within the next 10-25 years. 
The reason for that is the the orbital plane KOI-142b is inclined relative to the more massive 
KOI-142c, must precess around KOI-142c's orbital plane. As it precesses, 4 relative to the 
transit plane will increase to the point when bf, > 1. The exact time interval during which 
the transits of KOI-142b will still occur depends on the present v and Qc- The transit 
window can be as long as several decades for Qc — 90° and the smallest acceptable values of 
ic (~1.5°-2°). On the other hand, in can be as short as ~10 years if Qc — 270° and ic ~ 5°. 

The transit windows of KOI-142b should last 30-50 years, and should re-appear with 
the periodicity of Pq ~ 100 yr (Fig. |H]). It is unlikely that the transits of KOI-142c could 
ever be observed from the Earth, because the orbital plane of KOI- 142c remains inclined 
with respect to the transit plane (assuming that there are no additional massive bodies in 
the system). We find that 1.5 < 6c < 4 with the best fit values occurring for 6c — 3 (Fig. |H]). 

We predict that the RV measurements of KOI- 142 should reveal at least two basic 
periods. The RV term with a 22.34 d period, corresponding to KOI-142c, should have K ^ 
48 m s~^ amplitude (Table 4). The amplitude of the 10.95 d period term, corresponding to 
KOI-142b, is uncertain because we only have an upper limit on M^. With Mh/M^, = 5 x 10~^, 
we obtain i^ ~ 5 m s~^. 



4. Summary 

We performed a detailed analysis of the KOI-142 transit photometry. The observed 
TTVs with an amplitude over 10 hours are some of the largest ever detected for a KOI. 
We found evidence for periodic TDVs that are nearly in phase with the measured TTVs, 
but are much smaller in amplitude (~5 min). To our knowledge, this is the first time that 
TDVs were unequivocally detected for an exoplanet. A thorough dynamical analysis was 
then conducted to find plausible interpretations for the measured TTVs and TDVs. 

Based on this analysis, we confirm KOI-142 as a system of two planets just wide of the 
2:1 resonance (average Pc/Pb = 2.03). The inner Neptune-class planet, KOI-142b, shows 
transits at the present epoch. We predict that these transits should disappear in the next 
few decades. Monitoring the impact parameter (or equivalently the transit duration) of 
KOI-142b will help to better constrain the inclination and nodal longitude of KOI-142c. 
The outer sub- Jovian mass planet, KOI-142c, is currently not transiting, and we find that 
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the viewing geometry and dynamics preclude the visibility of KOI-142c transits from Earth 
in the future. 

The dynamical configuration of KOI-142b's and KOI-142c's orbits near the 2:1 resonance 
can be used to constrain the formation and migration history of the two planets (Beauge 
et al., in preparation). The relatively large orbital eccentricity of KOI-142c (mean Cc = 
0.055) cannot be explained by present gravitational perturbations from KOI-142b, or by tidal 
migration (Lithwick & Wu 2012, Batygin & Morbidelli 2012). It was probably established 
early, possibly during KOI-142c's formation in the proto-planetary gas disk (e.g., Lega et al. 
2013). 

We thank the Kepler Science Team, especially the DAWG, for making the data used 
here available. DN acknowledges support from NSF AST-1008890. DMK is supported 
by the NASA Sagan fellowship. JH acknowledges support from NSF AST-1108686. GB 
acknowledges support from NASA grant NNX09AB29G. 
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Table 1: Mean transit parameters derived by taking the average of each posterior sample 
across the M = 100 joint posteriors derived for each transit epoch. 



Parameter 


Value 


Pb [days 
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24.8+}i 


k 


n Q7O+0.028 

U.O(Z_Qg29 


^bU 


8q 1 40+0-085 


p*,6(circ) [gcm^^ 


9 An+0.41 


Tfe [hours 


o iq+0.19 

'-'•-'-'-'-0.16 


7i4,b [hours 


o 00+0.20 
3.2«_Qi7 


T23,6 [hours 


2.99+°;}^ 


Ti2,b - ^34,6 [hours 


0.142+°:°}! 



Table 2: Final values of the free parameters used in the dynamical fits of the transit times and 
durations (TTVs and TDVs). The upper limit of Mb/M^, corresponds to the 95% quantile 
of the posterior distribution. The orbital inclination i and nodal longitude Q are given with 
respect to the transit reference system (Nesvorny et al. 2012). 



Parameter 


Planet b 


Planet c 


Mp/M, 


< 5.2 X 10-5 


(6.321°:}^) X 10-^ 


Pp [days 

i[°] 

e 


10.954204+°-°™ 

0.931+°:°?^ 

0.05608l°:S 

90.9+?:I 


99 QQoq+0.0036 

ZZ.OOOO„QQQ3;^ 

3.8+11 

0.0564+°;°°^^ 

271.1111 


n[°] 


270 


25911^4 


A[°] 


R 419+0-094 


253.012-1 
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Table 3: Derived stellar parameters found by matching isochrones to the SPC-derived stellar 
effective temperature and metalicity as well as the light curve derived stellar density (ac- 
counting for orbital eccentricity) . For comparison, we show the Kepler Input Catalog (KIC) 
parameters derived by Brown et al. (2012). 



Parameter 


KIC 


This work 


Fitted param. 


M, [M©; 


0.96 


1 n99+0-023 

l.UZZ_QQ2g 


R* [Rq, 


0.74 


q61+°°2° 


TeS [K] 


5361 ± 135 


5513 ±67 


Fe/H] [dex 


-0.10 ±0.20 


±0.37 ±0.08 


log g [dex 


4.68 ±0.25 


4 489+0.018 


L* [Rq^ 


- 


n 7c;c:+0.056 

U.(00_QQg7 


Age [Gyr; 


- 


9 4c;+i-20 


Mv 


- 


tr ioq+0.116 
O.i00_o.o90 


Distance [pc 


- 


385_2g 



Table 4: Physical planetary parameters derived by combining the stellar parameters shown 
in Table and the dynamical parameters shown in Table [E The upper limits correspond to 
the 95% quantile of the posterior distribution. The orbital inclination i and nodal longitude 
Q are given with respect to the transit reference system (Nesvorny et al. 2012). 



Parameter 


Planet b 


Planet c 


Mp [Me] 


< 17.6 


215.9+^1 


Rp -Rffi 


4.23l°J° 


- 


pP [gcm-3; 


0.48+°i^ 


- 


Pp [days 


10.954204+°-°™ 


99 QQCQ+0.0036 
ZZ.OOOO_o.o031 


i[°' 


0.9621°:°?! 


3.71II 


e 


05611+°-°°i2i 


0.05628+°:°°^^ 


w{°' 


90.71?:^ 


271.1+1! 


n[°] 


270 


259+1^4 


A[°] 


fi 49c;+0.098 


253.0+^i 


i^[ms-i] 


< 5.0 


48.2lJ;i 
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Fig. 1. — A so-called "river-plot" for KOI-142.01, as used by Carter et al. (2012). Each row 
corresponds to a single transit epoch time series and the squares within that row represent 
the chronological sequence of flux measurements represented by the color. We denote a 
normalized flux of unity as gray and a normalized flux equal to a transit depth as black. 
White squares denote a lack of data. 
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Fig. 2. — A triangle plot for the dynamical parameters fitted to the transit times and 
durations of KOI-142b. Here, the variance covariances between each fitted parameter are 
easily seen. 
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Fig. 3. — A triangle plot for the dynamical parameters fitted to the transit times and 
durations of KOI-142c. Here, the variance covariances between each fitted parameter are 
easily seen. 
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Fig. 4. — The TTVs and TDVs of KOI-142.01. Data with associated error bars are shown 
in both panels. The TTV errors are ~l-2 min and are unresolved on the scale of panel (a). 
The green line shows our best simultaneous fit to the TTVs and TDVs. The red line shows 
the best TTV-only fit. These two fits give practically identical TTVs (the two lines overlap 
in panel a). 
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Fig. 5. — The short-period chopping of KOI-142b's TTVs produced by the orbital conjunc- 
tions with KOI-142c. The main resonant terms were removed from the measured and model 
TTVs by a Fourier filter. The red and green lines show our best TTV-only and simultaneous 
fits from Fig. HI 
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Fig. 6. — Orbital evolution of KOI-142b and KOI-142c. The red lines denote the orbital 
elements of planet b, the blue lines denote the orbital elements of planet c. The gray area 
in panel (d) shows an approximate region where transits occur for Q ^ 90° or i7 ^ 270°. 
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Fig. 7. — The predicted future behavior of TTVs (panel a), transit duration (panel b), and 
impact parameter (panel c) of KOI- 142b obtained from our best-fit solution. The two lines 
show a representative best solutions for i7 ~ 90° (green line, 4 = 2.7° at the reference epoch 
t = 0) and fi :ii 270° (red line, ife = 4.5° at t = 0). The TTVs in panel (a) were offset by 
±1000 min for clarity. The gray area in panel (c) shows the region where transits occur 
(6<1). 
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Fig. 8. — The same as Fig. [7] but for an extended time interval. The blue line in panel 
(c) shows the impact parameter of KOI- 142c corresponding to the best fit solution shown in 
Fig. EH. 



